The origin of normal heat conduction in one-dimensional classics system 
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We propose a new one-dimensional lattice model with strong asymmetric interaction potential 
and investigate heat conduction in this model numerically. We find that Fourier law is obeyed. 
Based on the phonon theory, we find a new scattering mechanism of phonon because of the breaking 
of the lattice segment. It is shown that in most of scattering process in this model momentum is 
destroyed as well as the Umklapp phonon-phonon scattering process which leads to the normal heat 
conduction. At last, we extend our analysis to the same class model with asymmetry interaction 
potential and get a general conclusion. 
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The heat conduction in one-dimensional classic sys- 
tem is a long debated problem in statistical physics. 
Although Peierls's theory [lj shows that the Umklapp 
phonon scattering process can lead to the normal heat 
conduction in lattice system, but it is not definitely clear 
in one-dimensional system what is the sufficient and nec- 
essary condition for the valid of Fourier law. This prob- 
lem has been attracting much attention in these two 
decades Q. It is previously believed that the break 
of momentum conservation may be the reason for the 
normal heat conduction. The numerical studies 
showed that in one-dimensional system subjected to on- 
site potential heat conductivity is independent on sys- 
tem size, which indicates the valid of Fourier law. The 
momentum is not conserved in these models because of 
external potential. However, such external potential is 
not nature in real material, just in some treatment of 
many-body system the mean field theory can produce 
this potential. On the other hand, it has been sug- 
gested analytically that the conservation of momentum 
is the key for the anomalous conductivity in low dimen- 
sion [6]. One can consider one-dimensional model with 
the mentum conservation without external potential such 
as Fermi — Pasta— Ulam(FPU) model. The interaction 
between phonons will arise because of nonlinear term in 
interaction between particles and it was expected there 
should be a finite heat conductivity. However, the nu- 
merical results showed that the heat conductivity k is 
size-dependent 7] and divergent versus the system size as 
k ~ L a , which was also confirmed by the mode-coupling 
theory kinetic theory Q and fluid hydrodynamic 
theory [l0| . The studies on one-dimensional hard-point 
model also showed the divergent of heat conductivity and 
it is believed that in this momentum conserved model 
heat conduction is anomalous and a is common univer- 



sity of this kind of model |ll|,|l2j. However, there are few 
results which are not in agreement with this belief. One- 
dimensional rotor model was showed to have a finite heat 
conductivity pjl EH • Recently, one-dimensional lattice 
system with asymmetric interaction potential was stud- 
ied numerically and found to have a finite heat conduc- 



tivity [15]. 

In this paper, our goal is to find the origin of normal 
heat conduction in one-dimensional system without in- 
troducing a external potential, We investigate the heat 
conduction in a new proposed one-dimensional lattice 
model using molecule simulation and find Fourier law is 
obeyed. Using the phonon theory, we find a new scat- 
tering process which does not occurs in one-dimensional 
lattice with symmetric interaction. It will destroy the 
momentum conservation of phonon and leads to the nor- 
mal heat conduction. 

We consider a new one-dimensional lattice system with 
Hamiltonian 
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where xi and pi are the position and momentum of ith 
particle respectively and V is the interaction potential 
energy, a is the lattice constant which is the equilibrium 
distance between pairs particles. We take a special form 
of the interaction potential. We assume the interaction 
between pairs particles is still spring like. However, when 
the distance between neighbor particles is smaller than 
the lattice constant, there is interaction between them. 
If the distance is larger than lattice constant, the inter- 
action vanish. Hence, this interaction potential is 
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In Fig.l, we give a schematic picture of the interaction 
potential. We can find this model has not only the prop- 
erties of lattice, but also of the gas model. 

To investigate the heat conduction using the computer 
simulation, the model should be connected with heat 
bath. We adapt Noose-Hoover [l(| heat bath on the ends 
of the present model. By setting the temperature differ- 
ence, we can calculate the heat current in this model after 
the thermal steady state is approached. 

We find that in this model temperature gradient can 
be well formed along the chain if the simulation time is 



2 




-1 1 500 1000 1500 , 2000 2500 3000 



FIG. 1: The schematic of the interaction potential energy. 



long enough, which indicates the local thermal equilib- 
rium state can be reached. In Fig. 2, we plot the heat 
conductivity K as the function of the system size L, we 
note that when the system size is small, the heat con- 
ductivity is divergent. However, when the system size is 
larger than 2000, k become divergent slowly and comes 
into a saturate value which mean we can expect a nor- 
mal heat conduction. This result shows that even in a 
momentum conservation one-dimensional system, a hnitc 
heat conduction can be found. 

Comparing the interaction potential used here with the 
interaction potential used in ref . , we can see although 
they have different form, but they both have asymmetric 
properties. And, the asymmetry in our model is much 
stronger than others and heat conductivity in our model 
is easy to converge. So a question arises, is the asymmet- 
ric of interaction lead to normal heat conduction? If so, 
what is the underlying mechanism. In the following, we 
will investigate these problems. 

First we consider one-dimensional harmonic chain of 
particle each of them has mass m and is coupled each 
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other by the interaction V(x) = -j-x 2 . With the period 
boundary condition a set of normal mode coordinates are 
given by 



Qk = —j= ^2 x i ex P{ikl), 
1 N 

Pk = —f=^2piexp(ikl), (3) 
v N i=i 

where k is the wave vector and given by k — 2nj/N 
and j is a integer chosen so that k is in the Brillouin 
zone (— 7r, 7r). In quantum mechanic, one can use the 
annihilation and creation operators at and a^~,and a k is 



FIG. 2: The heat conductivity k as the function the system 
size. The temperature at two heat bath are T+ = 0.3 and 
T_ =0.2 respectively 



defined as 

where flk are the normal mode frequencies given by 

k 

fl k = 2uj | sin - \=uj a uj k , (5) 
the Hamiltonian is diagonalized as 

k 

The quantum of vibrational energy is called a phonon. In 
the Peierls theory of solid, the heat conduction is carried 
by phonons, if there is no scattering between phonons, 
the mean free path will become infinite and temperature 
gradient can't formed. By introducing the nonlinear in- 
teraction between particles, it is shown that there exists 
collisions between phonons. If the momentum of phonons 
is conserved during the collision, this process is called 
normal process and if the momentum is not conserved it 
is called Umklapp process. Only the later will contribute 
to the normal heat conduction because in which the heat 
resistance can be formed. In the previous works on this 
subject, the main aim was to find a scattering mechanism 
to cause normal heat conduction. The external potential 
is believed to be a key to answer this question. In real 
space the momentum conservation is destroyed by the 
external potential, so it is suggested that it is true for 
phonon. 

Let's consider our model which has the features of both 
the lattice model and gas model. We can not apply 
phonon theory to it directly because pair particles have 
no interaction if the distance between them is larger than 
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ftfc^. In the breaking of segment, it will split into many 
phonons. The wave vectors of new generated phonons 
take the values in equation (JSJ). Consider a simple three 
phonons process that one phonon splits into two phonons, 
we denote their momenta as hk^ and ftfc 3 ", it is easy to 
find that the momentum is not conserved that 



FIG. 3: The schematic of the breaking of lattice segment 



a. Assuming the system is composed of N particles and 
the length is L — Na with the period boundary condi- 
tion, when it is in ground state that each particle is at 
it's equilibrium position, there is no force between them. 
When an excitation emerges along the chain, particles are 
pushed from their equilibrium position. Some of them 
are close to its neighbor and the interaction occurs be- 
tween them. As a result, a lattice segment is formed and 
between segments the distance is larger than a, which 
means the model in thermal state is composed by lattice 
segments. We can use phonon theory to describe lattice 
segment. 

Assuming a lattice segment composed by N± particles 
and its length is L\ — N\a. The wave vectors of phonons 
excited in it are 

2tt 

ki = -j-ji, (7) 

where j is a integer and takes the values from 1 to N% — 1. 
In the subsequent evolution of model, it is expected that 
at certain time, there will be a breaking in this segment 
that one particle will move far away one of its neighbor 
and their distance is larger than a. As a result, the orig- 
inal segment now splints into two segments as shown in 
Fig. 3. Assuming the number of particles in these two 
segment are N2 and N$ and the lengths of them are L2 
and L3 respectively. Hence, the wave vectors of phonons 
in these two new lattice segments are 

k 2 = ~rj2, 

2tt 

k 3 = ~rJ3, (8) 

where j% and j'3 are also positive integers and their max- 
imal values are N2 — 1 and A3 — 1. As L\ = L2 + L3, it 
is easy to see 

k 2 > h,k 3 > k%. (9) 

After the breaking of segment phonons are changed. 
That is to say that the breaking of segment is a mech- 
anism that causes the change, namely scattering, of 
phonon. Let's consider the one in original segment with 
the wave vector k[ = 2n/Li and its quasi- momenta is 



fik[ ^ hk'i + Kk'l', (10) 

due to the equation®, k' 2 ' and fc 3 " are larger than k[. To 
make the equivalent of both sides of above equation, one 
should add another momenta into left side. We denote 
its wave vector by G and get a equation as following 

hk[ + HG = hk'2 + hk'z'. (11) 

This is a Umklapp process. In fact, in the breaking of seg- 
ment, the phonon k' x should split into many new phonons 
in two new segments. The wave vector of each new 
phonon is larger than the old one. So, it is also showed 
that such multi-phonon process should satisfy equation 
as 

mi 7722 

hk^ + hG' = J2 hk 2 + J2 hk 3> ( 12 ) 

i=l 3=1 

where mi and 7712 are the numbers of phonons involved in 
this process in the two new segments. G' is an additional 
term to fulfill the equality of equation. This process is 
also a kind of Umklapp process. For other phonons in 
initial segment, it can also be showed in this way that 
in most scattering processes the momentum conserva- 
tion are destroyed, although for the phonons with larger 
wave vector in segment before breaking, there will be few 
processes that satisfies momentum conservation. On the 
other hand, as is well known, only the phonon with low 
wave vector contribute to the heat conduction. In term of 
our illustration, the scattering processes of phonons with 
small wave vectors are almost Umklapp process. Hence, 
we can conclude that the breaking of lattice segment in 
this model is a scattering mechanism of phonons which 
destroys the momentum conservation. As a result we can 
expect the normal heat conduction in this model con- 
firmed by the numerical result in Fig. 2. 

As our model has the asymmetric interaction which 
cause the breaking in the evolution of model, next, we 
extend our analysis to the models with other asymmetric 
interaction, such as Lennard-Jones potential written as 

V(r)=e[C—) 12 -2C—)% (13) 
r r 

where r is the distance between pairs particles. For one- 
dimensional lattice, r — Xi — Xi~\. r m is the distance at 
which potential has the minimum, e is the depth of po- 
tential well. As well known, Lennard-Jones potential is 
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are defined as following 




FIG. 4: The rate 7 versus the relatively difference of displace- 
ment from the minimum Ar in the Lennard- Jones potential, 

£ = 1.0. 



asymmetric around the minimum. To show the asymmet- 
ric more clearly, we calculate the differential of potential 
around the the minimum of potential and define a rate 
as 



dV{r)/dr \ r =r m +Ar 
dV(r)/dr \ r =r m -Ar 



(14) 



where Ar is the relative displacement of particle from 
its equilibrium position. Some results are illustrated in 
Fig. 4, we can find that for fixed parameters of poten- 
tial the differential of potential at the same displace- 
ment on the two sides of minimum become more different 
when the displacement Ar becomes large, which shows 
the asymmetric of interaction potential. And, we find 
that r m determines the degree of asymmetry of poten- 
tial. When r m becomes small, the differential of poten- 
tial, namely, the force acted on pairs particles becomes 
much smaller when the two particles leave apart than the 
case when two particles move close to each other. 

To use phonon theory, the general treatment of nonlin- 
ear potential is to approximate it by expanding it around 
the minimum and just take the second order term as the 
first order term is zero. However, as the expansion co- 
efficient of second term, dV/dr is taken at the minimum 
no matter what sigh of Sr. Hence, the asymmetry is lost 
in the approximation. 

To involve the asymmetry in the linear approximation, 
we expand the potential along the trajectory. We can get 
the function of linear term. For simplicity, we use the 
average linear term function. And the form of it on the 
left side and right side of minimum denoted by ki and k r 
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Considering the case r m — n — r r — r m , it is easy to 
conclude from Fig. 4 that ki > k r , showing the asymme- 
try. Thus, we can consider the evolution of this model 
in linear approximation. The coupling constant between 
pairs particles are dependent on distance, taking value 
ki or k r . In other words, it is time-dependent. In the 
evolution, with the change of distance of pairs particles, 
there will be a crossover at which the couple constant 
will switch between fc; and fc r , then the phonon spec- 
trum will change. It is a scattering mechanic of phonons. 
Consider the small r m , it can be expected that ki 3> k r , 
we can ignore k r and only ki exists, then the model with 
Lennard-Jones potential is reduced to our model. Hence 
our finding of the scatting of phonon can be applied to 
such model and the same conclusion is obtained. 

In conclusion, we found normal heat conduction in our 
one-dimensional lattice model with asymmetric interac- 
tion, that provides another example which obey Fourier 
law with momentum conserved. It is to say that the 
momentum conservation is not the necessary condition 
for normal heat conduction. Furthermore, we found in 
our model, due to the asymmetric interaction potential, 
there is a scattering mechanism caused by the breaking of 
the segment, which leads to the normal heat conduction. 
That is to say the phonon spectrum is time-dependent. 
We can apply our analysis to other model with asymmet- 
ric interaction. As is well known, the asymmetric inter- 
action potential is the reason for the expansion of solid 
potential 17]. Our finding indicates that it can account 
for the normal heat conduction. 
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